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Abstract 



For the vertex corrections to the partial decay rate T(Z — > bb) involving the top 
Q_i! quark only the leading terms of order aa s in the 1/-M* expansion are known. In 

this work we compute the missing next-to- leading corrections. Thus at 0(aa s ) the 
complete corrections to the decay of the Z boson into bottom quarks are at hand. 



PACS numbers: 12.38.Bx, 13.30.Eg, 13.38.Dg, 14.65.Fy 



1 Introduction and notation 

At the Large Electron Positron collider (LEP) at CERN approximately four million decays 
of the Z boson per experiment have been observed. Because of this enormous statistic 
a lot of observables have been measured with very high precision - - sometimes in the 
region of a few per mille or even below. Also the hadronic decay of the Z boson has been 
measured to an accuracy of roughly 0.1% and amounts to Thad = 1742.4 ± 3.1 MeV ]I[. 

From the theoretical side also a lot of effort has been undertaken in order to give a 
precise prediction for Thad- The QCD corrections are known up to 0{a z s ) both in the 
massless limit g] and at 0(M 2 /s) || §J. A complete 0(a 2 s ) result [|] and the leading 
terms at 0(a^) |j] are available for the singlet contribution (for reviews see also |7|, ||). 
Concerning purely electroweak corrections a complete calculation is available at one-loop 
level 0, whereas at two loops the leading terms in M t are known [Kj. The mixed 



corrections of 0(aa s ) for the decay of the Z boson into the quark flavours u, d, s and c 



were considered in [11], [12], |13| . From these results also those 0{aa s ) contributions to 
YiZ — > bb) can be extracted where a photon or a Z boson is exchanged between the 
bottom quarks. For the other class of diagrams contributing to the decay into bottom 
quarks, namely those involving W bosons and top quarks in the loop, only the leading 
M 2 corrections [14, 15 1 and the hiM 2 terms [16, 17] are at hand. In this work we close 



the gap and provide the constant term at next-to-leading order. In addition three more 
terms in the high-M 4 expansion are computed. It will be demonstrated that these five 
terms provide a reliable result with negligible errors. 

It is useful to distinguish in the decay of the Z boson between universal and non- 
universal corrections. Universal terms are independent of the produced fermions and 
arise from corrections to the gauge boson propagators. The non-universal corrections, 
sometimes just referred to as vertex corrections, depend on the fermion species considered. 
Of special interest is thereby the decay into bottom quarks since in this case an additional 
dependence on the top quark mass appears. 

The partial decay rate of the Z boson into bottom quarks can be written in the form 

r(z-*bb) = r°{v 2 b +al)+5r univ + 5r b , (i) 

with r° = N c M z a /12s|c|, v b = —1/2 + 2s|/3 and a b = —1/2. s e is the sine of the weak 
mixing angle and c@ = 1 — s@. <5r un i v represents the universal corrections and contains 
contributions from the transversal part of the renormalized Z boson polarization func- 
tion evaluated for q 2 = M§, Ar, the radiative corrections entering the relation between 
Gp,C(,Mz and M^y, and the universal corrections arising from the 7 — Z interference. 
5T b comprises all terms directly connected to the Zbb vertex. This includes both pure 
QED, QCD and electroweak corrections and the mixed contributions of order aa s . In 
this work we will concentrate on non-universal corrections to ST b where the top quark is 
involved, in the following denoted byQ ST^ . The corresponding diagrams are pictured in 

1 Although the index W has been chosen, in a covariant gauge also the diagrams involving the charged 
Goldstone bosons 6^ have to be considered. 



Fig. 0. To get the 0(aa s ) corrections an internal gluon has to be attached to them in all 
possible ways. Analogously, we define 5T b as the contribution arising from the diagrams 
with internal Z boson exchange. The mass of the 6 quark will be neglected throughout 
this paper. 

Whereas ST b is finite after the Z boson contribution of the wave function renormal- 
ization of the quarks is taken into account, for 5T^ an additional counterterm induced 
by the Born and (D(a s ) result has to be added. In Feynman gauge it reads^ jl2| : 



5Ti 



srf - 


= dT b +dT b , 


^ct,W 


Sq it °U 


1 1 2 , a s n ( 1 

.-6-3 ce+ v c n-8 




5 5 2 / 1 2 ,\ , / 1 

-18-9 Ce+ v-3-3 Ce r^ + l-6 




71 


55 55 9 / „ ,\ A 
.-48-24 c|+ ( 1 + 24 ) 6 




+ {-\-\ Cl ) l » W+ {-\-\ Cl ) U 


320569C 


3, l,M = ln(// 2 /M^) and l e = mc|. 5T° b ' w 



-e 



e 



'e 



e 



(2) 



with £ 3 ps 1.202056903, /^vk : = ^(/i^/M^) and /e = lnc@. 5T b ' yy is the contribution from 
the Feynman diagrams depicted in Fig. ||. 

In principle there are several options to write down an equation for T(Z — > 66) con- 
taining the radiative corrections. Eq. (||) is only one possibility. Often parts of the vertex 
corrections are comprised into so-called effective couplings and the QCD corrections are 
described by a factor (1 + a s /n). This will be discussed in Section | where the results are 
presented. In the next Section some details on the calculation are given. 



2 The calculation 



In the following the method used for the calculation is briefly described. Instead of 
directly computing the two-loop vertex diagrams the three-loop polarization function of 
the Z boson is considered. The imaginary part then immediately leads to contributions 
to ST^ . This avoids the separate treatment of infra-red and collinear divergences and has 
the advantage that the advanced computational tools available for two-point functions can 
be used. In addition, as the top quark is much heavier than all other mass scales involved 
in the problem, it is tempting to perform an expansion in the inverse top quark mass. 
The method which provides a systematic expansion in 1/M t is given by the so-called hard 
mass procedure (HMP) [|T8[|. 

The application of the HMP requires in a first step the identification of all hard sub- 
graphs which contain all lines carrying the large mass and whose connectivity components 
are one-particle-irreducible with respect to lines corresponding to light or massless par- 
ticles. These subgraphs have to be expanded with respect to the small masses and the 



2 We are using dimensional regularization with space-time dimension D 
the ln47r and je terms in all formulas. 
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Figure 1: Diagrams contributing to 5T^ in 0(a). Thin lines correspond to bottom quarks, 
thick lines to top quarks, dotted lines to Goldstone bosons and inner wavy lines represent 
W bosons. 



external momenta. In the full diagram all lines of the hard subgraph have to be shrunk 
to a point to give the so-called co-subgraph. The result of the expansion is inserted as an 
effective vertex. Finally the loop integrations are performed. 

At two-loop level seven diagrams have to be considered. Their calculation would still be 
feasible by hand. At three loops, however, 69 diagrams contribute and a computation by 
hand is very painful especially if higher order terms in the 1/M t expansion are considered. 
The HMP applied to the 69 initial diagrams contributing to 5T^ results in 234 sub- 
diagrams which have to be expanded in their small quantities. For this reason the program 
package EXP written in Fortran 90 was developed. It generates all possible subgraphs 
of a given diagram together with some auxiliary files which allow the computation to be 
done automatically. The output of EXP can directly be used as input for MATAD jjHJ 
and MINCER [f2(J where the loop integration is performed. 

It is possible to separate the diagrams at two-loop level (see Fig. ||) into two classes: 
If the HMP is applied to the diagrams of type (a), (b), (d) or (e) the resulting integrals 
factorize into massive diagrams with only one mass scale times purely massless ones 
with external momentum q. For the types (c), (/) and (g) the HMP also leads to a 
factorization of the original integral. Here, however, massive one-loop integrals with two 
different masses and non-vanishing external momentum have to be computed. At three- 
loop level, where the same classification is valid, in principle even the 0(e) part would 
be necessary. For this reason we apply the HMP again to these co-subgraphs with the 
conditions M^ » iw^-w ^ 9 2 ' or alternatively, iw^-w ^ ^w ^ Q' '■ Here, £ w is the 
gauge parameter appearing in the W and 4> propagators. Of course, both descriptions 
must lead to identical final results as the £w dependence drops out at the very end. In 
intermediate steps, however, the expressions which have to be evaluated are different. 

At the end q 2 = M§ has to be chosen. Then the above inequalities are seemingly 
inadequate. However, a closer look to the corresponding diagrams shows that actually 
(2M W ) 2 , respectively (M w + M t ) 2 , has to be compared with q 2 . The procedure is further- 



more justified a posteriori as rapidly converging series are obtained from the expansions. 

3 Results 

Since we are interested in the virtual effect of the top quark which renders the decay of 
the Z boson into bottom quarks different from the one into other down-type quarks, in 
the following the quantity 5T b ' — ST^ is considered. In this difference the counterterm 
contribution exactly cancels. This means that the renormalized partial decay rate can be 
written in the form 



sr? = srY+(srf w -s^ 



w 



(3) 



where 5T^ is the contribution from the diagrams involving a W boson to the partial decay 
rate T(Z — > dd) [fL2| . Note that besides the 1/e poles also the dependence on £ w drops 
out in the difference ST k ' — 5Tj . For convenience the following notation is introduced: 
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For STl_ d we get: 
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with y = 1/(44), C2 = tt 2 /6, L tw = ln(M?/M w ) and l @ as in Eq. (@). Note that the 
successive application of the HMP to the diagrams of class (c), (/) and (g) (see Fig. |2|) 
leads to an expansion in 1/cq. However, the coefficients decrease very rapidly which 
justifies this strategy. Moreover, a W boson in the final state is always accompanied by 
either another W boson or a top quark. For this reason we choose the variable y for 
the presentation of the results. Nevertheless, two remarks concerning the expansion in 
1/cq are in order: First, note that the above series is nested and therefore each order in 
1/M t 2 produces an additional power in y. To demonstrate the quality of the convergence 
all available terms are displayed. Second, one recognizes that the coefficients of the 
logarithms seem to be truncated series in l/c@, which suggests that they might be exact. 
These remarks are also taken over to the 0(aa s ) corrections given by 
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For the definition of the top quark mass, M t , the on-shell scheme has been adopted. With 
the help of the equation 



M t 



m t {lA 



Ti \ 4 mt(fi) : 



(7) 



the transformation to the MS scheme can be performed. 

There are several checks for the correctness of our results. First of all, our new calcu- 
lation with the automated HMP provides an independent check of the M| and the In M 4 2 
term. The terms up to (and including) order 1/M t 2 and 1/cq were calculated with arbi- 
trary gauge parameters £w and £5 for the electroweak sector and QCD, respectively. We 
could check that to each order in the 1/Mt expansion £w drops out separately. Whereas 
for the M t 2 and 1/Mf terms this happens after taking the sum of all contributing diagrams, 
for the (M 2 ) order only the difference 5T b ' — 5T d ' renders the result gauge invariant^. 
£s already drops out if only the sum of diagrams contributing to the separate classes (see 
Fig. HD is considered. The 1/Mf and 1/Mf corrections and the higher order terms in 
the 1/cq expansion were computed in Feynman gauge as then the calculation is much 
faster. Note, that in Eq. (||) where all parameters are expressed in the on-shell scheme 
the explicit \x dependence drops out which is also a welcome check for our calculation. 

3 It should be noted that we partly had to repeat the calculation done in |12] for arbitrary gauge 
parameter £w 



Let us now discuss the numerical relevance of the newly computed terms. Using 



Sq = 0.223 the corrections induced by PU-bosons read in the on-shell scheme: 
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Using Eq. (0) one may transform the top mass, M t , to the MS scheme. For the renormal- 
ization scale /i explicitly appearing in the resulting expression, 5T ^ , we adopt the choice 
fi 2 = M z . One finds {m t denotes the MS top mass, and l t w = m ( m ?/^vlO) : 
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One observes that for realistic values of Mz and M t , respectively m t , the constant at 
next-to-leading order dominates over the InM 2 term known before. Furthermore, the size 
of the contributions from the individual 1/M 2 orders is roughly the same in both schemes. 
It is remarkable that at one- loop level the corrections arising from the 1/M 2 terms are 
of similar size than the one from next-to-leading order, however, the signs are different. 



The higher order corrections in 1/Mt are smaller, which means that effectively only the 
leading M t 2 term remains. Proceeding to two loops the situation is similar: Starting at 
0(1/ M t 2 ) the sign is opposite as compared to the leading terms and a large cancellation 
takes place. Here, the 1/M t term is still comparable with the 1/M 2 contribution. The 
1/M f 6 term, however, is small and thus strongly suggests that the presented terms should 
provide a very good approximation to the full result. 

For M t = 175 GeV, M z = 91.19 GeV, a = 1/129 and a s (M z ) = 0.120 we get: 

8TY_ d = (-5.69-0.79 + 0.50 + 0.06) MeV = -5.92 MeV. (10) 

The first two numbers in Eq. fllPf) correspond to the 0(a), the second two to the 0(aa s ) 
corrections. Each of these contributions is again separated into the M 2 terms and the 
sum of the subleading ones. For the diagrams containing a Z boson |T^| the numerical 
value reads: 

(JTf = (0.52-0.01) MeV = 0.51 MeV. (11) 

One can see that even the contributions of the subleading terms of diagrams arising from 
a W boson exchange are more important than the corrections resulting from 8T% . 

Let us finally also give the net effect of the non-factorizable corrections. Adopting a 
notation analogue to Eq. (U) and using the expansion for ST^ in the limit of small y [|12|Q 
the difference to the result where "naive" factorization is assumed reads 



^(C F 5rP' w -8T^ w ) = 0.68 MeV. (12) 

This is comparable both to the error on the present experimental value of r^ and even 
to the one-loop corrections from the Z boson diagrams (see Eq. (|TT|)). For T(Z — ► bb) the 
factorization is often performed such that the leading M 2 term is reproduced correctly. If 
this is taken into account the deviation corresponding to Eq. ( |P2"|) reduces to —0.04 MeV 
which is still larger than the mixed 0(aa s ) corrections induced by internal Z bosons 
(second number in Eq. (|TT|)). 

To conclude, the missing non-universal piece to the decay of the Z boson into bottom 
quarks to 0(aa s ) has been computed. An expansion for large top quark mass has been 
performed and it was demonstrated that only the inclusion of power suppressed terms in 
1/M 2 leads to reliable predictions. The results are presented in a form which allows a 
simple implementation into program libraries for the description of the Z line shape (see 
and references therein). 
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